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ABSTRACT 


Three  additional  elements  have  been  added  to 
the  element  library  of  the  FINEL  finite  element 
program:  a rectangular  membrane  element,  a rect- 

angular combined  membrane  and  bending  element,  and 
a triangular  combined  membrane  and  bending  element. 

This  report  describes  the  characteristics  of  the 
three  new  elements  and  includes  data  on  the  bending 
convergence  behavior  of  the  combined  elements.  The 
stiffness  matrices  of  the  three  elements  are  pre- 
sented in  the  appendix. 

ADMINISTRATIVE  INFORMATION 

This  project  was  conducted  at  the  Naval  Ship  Research  and 
Development  Center  as  part  of  a study  on  computer  methods 
applied  to  naval  ship  structures.  Funding  was  provided 
through  NAVSHIPS  under  Subproject  SF  35  422  302,  Task  11989. 

INTRODUCTION 

In  view  of  the  popularity  of  the  FINEL  finite  element 
1 2 

program,  ’ an  effort  was  made  to  increase  its  versatility  and 
efficiency  by  adding  several  elements  to  the  finite  element 


References  are  listed  on  page  39  . 
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library.  Three  types  of  elements  were  chosen:  a rectangular 

membrane  element,  a rectangular  combined  membrane  and  bending 
element,  and  a triangular  combined  membrane  and  bending  ele- 
ment. The  combined  membrane  and  bending  elements  were 
selected  because  the  original  FINEL  program  contained  no  plate 
bending  elements  and  thus  could  not  be  applied  to  a large 
class  of  problems  that  are  both  interesting  and  significant. 
The  rectangular  membrane  element  was  added  to  the  element 
library  because  of  its  potential  for  increasing  the  operating 
efficiency  of  the  program.  The  original  program  already  in- 
cluded a general  quadrilateral  membrane  element,  but  past 
experience  with  FINEL  had  indicated  that  this  quadrilateral 
element  was  most  frequently  used  with  rectangular  geometry. 
Thus,  the  internal  mathematics  involved  justified  inclusion 
of  a rectangular  membrane  element  as  a separate  element  in  the 
library . 

The  most  important  characteristics  of  the  specific  three 
new  elements  are  described  herein. 

STIFFNESS  MATRIX  DEFINITIONS 

Before  describing  the  specific  elements,  it  is  necessary 
to  discuss  briefly  two  somewhat  different,  but  entirely  con- 
sistent, definitions  of  the  element  stiffness  matrix  because 
reference  will  be  made  later  to  both  of  these  definitions.  An 
understanding  of  basic  finite  element  theory  is  assumed  here; 
some  readers  new  to  this  technique  may  wish  to  first  consult 
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one  of  the  more  comprehensive  sources  available  on  the 


, 3,4 

theory. 

The  local  element  coordinate  system  and  the  sign  conven- 
tion employed  in  this  report  and  in  the  FINEL  program  are 
illustrated  in  Figure  1.  A right-handed  system  is  adopted 
here.  Translational  displacement  and  force  components  are 
positive  in  the  positive  direction  of  the  corresponding  ele- 
ment axes  whereas  rotational  displacement  and  moment  compo- 
nents are  positive  in  the  usual  right-handed  sense.  The 
coordinate  system  and  sign  convention  differ  in  certain  re- 
spects from  those  employed  in  various  references  cited  herein. 
In  view  of  all  the  ramifications  which  arise  because  of  these 
differences,  extreme  care  when  consulting  these  (or  any) 
references  on  finite  element  theory  cannot  be  overemphasized. 

Consider  some  general  finite  element  that  possesses  a 
total  of  m nodal  displacement  components  and  n rigid  body 
degrees  of  freedom.  Since  the  node  displacements,  in  general, 
consist  of  the  superposition  of  a rigid  body  displacement  of 
the  entire  element  and  the  deformation  of  the  element,  it  is 
convenient  to  divide  the  complete  column  matrix  of  node  dis- 
placements u into  two  smaller  column  matrices 


u 


(1) 


(2) 


(1) 


(2  ) 

where  u is  a matrix  of  n node  displacements  properly  chosen 
to  express  the  rigid  body  displacement  of  the  element,  and 
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u 1 ^ is  the  matrix  of  the  remaining  m-n  displacements.  The 
complete  set  of  node  forces  represented  by  the  matrix  R may  be 


similarly  divided  into  two  smaller  matrices  R ^ ^ and  R^^ 


R = 


(2) 


where  the  components  of  R ^ ^ and  R ^ ^ are  in  one-to-one  corre- 
spondence with  the  displacement  components  of  and 


respectively.  Reference  2 defines  an  element  de  f o rma  t i on  ma- 


trix v which  is  equal  to  the  difference  between  the  total  u 


(1) 


and  that  part  of  resulting  from  the  rigid  body  displace- 

( 2 ) 

ment  (represented  in  terms  of  the  matrix  u ).  Thus  the  de- 
formation matrix  v can  be  related  to  the  node  displacements  by 
means  of  a matrix  P such  that 

v = P 

= P • u 

Using  the  matrix  defining  the  element  deformation,  the 
stiffness  matrix  k is  defined  according  to  the  matrix  equation 


R 


(1) 


= fe 


v 


(4) 


This  first  definition  of  a stiffness  matrix  is  the  one  em- 
ployed by  the  FINEL  program.  This  is  not,  however,  the  defi- 
nition most  often  encountered  in  the  literature  on  finite 
element  theory.  Nevertheless,  as  will  be  demonstrated,  the 
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second  definition  can  be  derived  in  a straightforward  manner 
from  Equation  (4). 

It  can  be  easily  shown  that  the  nodal  force  matrices  R 

and  are  related  by  the  matrix  P according  to  the  expres- 

sion 

T ( 1 ) 

R = P • R1,  ' (5) 

where  the  superscript  T indicates  the  matrix  transpose.  If 
Equation  (3)  for  v is  substituted  into  Equation  (4)  and  both 

T 

sides  of  the  resulting  expression  are  premultiplied  by  P , the 
following  relationship  is  obtained: 

T ( 1 ) T 

P R^l;  = P * fe  • P • u (6a) 

Using  Equation  (5),  the  above  may  be  rewritten  as 

R = PT  • fe  • P • u (6b) 

By  combining  terms,  this  expression  may  be  further  rewritten 
as 

R = k • u (7) 

where  the  matrix  k is  defined  as 

k = PT  • fe  • P (8) 

It  is  now  easily  seen  that  the  matrix  k defined  by  Equa- 

★ 

tion  (7)  also  represents  a stiffness  matrix.  This  second 
definition  of  the  stiffness  matrix  is  the  one  which  appears 


* 

The  matrix  k defined  here  should  not  be  confused  with  the 

similar  symbol  used  in  Reference  2 since  in  this  case  their 
definitions  differ. 
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more  frequently  in  the  literature  on  finite  element  theory. 

The  simple  relationship  between  the  two  stiffness  matrices  k 
and  k is  given  by  Equation  (8).  This  expression  indicates  how 
the  matrix  k may  be  calculated  from  a knowledge  of  k (and, 
obviously,  P) , but  it  can  also  be  shown  that  fc  may  be  obtained 
from  k simply  by  inspection. 

THE  THREE  NEW  ELEMENTS 

THE  RECTANGULAR  MEMBRANE  ELEMENT 

The  rectangular  membrane  element  is  an  isotropic  element 
of  constant  thickness  t with  linearly  varying  normal  stress 
distributions  along  the  alement  axes  and  a constant  shear 
stress  distribution.  This  element  and  its  orientation  with 
respect  to  the  local  element  coordinate  system  is  shown  in 
Figure  2 . 

This  rectangular  element  represents  a special  case  of  the 
general  quadrilateral  membrane  element  included  in  the  library 
of  the  original  FINEL  program  (as  Element  2)  and  described  in 
Reference  2.  The  rectangular  geometry  results  in  certain 
mathematical  simplifications  which,  at  least  in  theory,  pro- 
duce an  element  whose  stiffness  matrix  can  be  calculated  with 
greater  speed  and  accuracy  than  that  of  the  general  quadri- 
lateral . 
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The  stiffness  matrices  for  both  the  general  quadrilateral 
and  the  rectangular  membrane  element  are  defined  in  terms  of 
node  displacement  and  node  force  matrices  given  by 


u = 

{uiviVjufcvfcuivi} 

(9a) 

u(1)  - 

{ujufevfeulvl} 

(9b) 

u(2)  = 

{ U . V , V . } 

1 i J 

(9c) 

and 

R = 

{XiYiV;jXfeWl} 

(10a) 

R(1)  = 

( X j Xfe  Yfe  X 1Y  i } 

(10b) 

r(2)  = 

{X.Y.Y. } 
i i 3 

(10c) 

Employing 

the  customary  energy  approach,  it  is 

2 

sh  own that  the 

st if  fness 

matrix  k is 

computed  from  the  matrix 

product 

k 

T 

- 1 - 1 
= (B  ) • J • B 

(ID 

where  the  matrices  B and  J are  functions  of  the  element  geome- 


try and  material  properties.  These  matrices  are  presented  in 
Reference  2 for  the  general  quadrilateral. 

For  the  general  quadrilateral  element,  the  form  of  ma- 
trix B (a  square  matrix  of  rank  5)  is  sufficiently  complex  to 
make  manual  inversion  as  required  by  Equation  (11)  impracti- 
cal. Consequently,  the  stiffness  matrix  for  this  element  is 
calculated  in  the  FINEL  program  by  numerically  inverting  the 
matrix  B and  performing  the  necessary  matrix  multiplications 
within  the  program. 


A 
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The  simplified  geometry  of  the  rectangular  membrane  ele- 
ment allows  for  more  direct  calculation  of  the  stiffness  ma- 


trix. The  matrix  B in  this  case  is  quite  suitable  for  manual 
inversion  and  thus  it  is  possible  to  perform  the  necessary 
matrix  multiplications  algebraically  to  determine  in  closed 
form  the  individual  components  of  the  stiffness  matrix. 

Matrices  B,  B ^ , J,  fe  , P,  and  k (in  symbolic  form)  for 
the  rectangular  membrane  element  are  shown  in  the  Appendix  to 
this  report. 

THE  RECTANGULAR  COMBINED  MEMBRANE  AND  BENDING  ELEMENT 

Like  the  pure  membrane  element,  the  rectangular  combined 
membrane  and  bending  element  is  an  isotropic  element  with 
uniform  thickness  t.  The  stiffness  matrix  for  this  element  is 
obtained  by  combining  the  stiffness  matrices  of  a rectangular 
membrane  element  and  a rectangular  bending  element.  Since 
this  indicates  that  the  membrane  and  the  bending  behavior  are 
assumed  to  be  uncoupled,  the  two  forms  of  behavior  may  be  con- 
sidered separately.  The  Appendix  describes  the  manner  in 
which  the  stiffness  matrices  corresponding  to  the  two  forms  of 
behavior  are  combined  to  form  the  combined  stiffness  matrix. 

The  rectangular  membrane  element  chosen  to  provide  the 
membrane  stiffness  matrix  is  the  element  discussed  in  the 
immediately  preceding  section.  Consequently  the  membrane 
characteristics  of  this  element  and  those  of  the  combined 
membrane  and  bending  element  are  identical.  Since  the 
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characteristics  of  this  membrane  element  have  already  been 
adequately  discussed,  the  remainder  of  this  section  will  be 
devoted  mainly  to  a description  of  the  rectangular  bending 
element  selected  to  provide  the  bending  stiffness  matrix. 

For  the  combined  membrane  and  bending  element,  the  stiff- 
ness matrices  k and  k are  defined  in  terms  of  the  node  dis- 
placement and  node  force  matrices  given  by 


u = {u1viw1a161ujv:!w:ja:)6:jufevfewfeafeBfeu1v1wiai81} 


(1) 


{u.w.a.8.u,v,w,ei,8(U,v,w,oi,  8,  } 

Jjjjfcfcfcfefc  11111 


u(2)  = {u.vwaB.v.} 
l i i i l j 


and 


f - <v1^^vjVjVfcVfe,w1v1z1K1v 

R(1)  . (xJV;lVfeVfe''kV1VlK1V 
R(2)  = {X.YiZiK1MiYj } 


(12a) 

(12b) 

(12c) 

(13a) 

(13b) 

(13c) 


In  these  expressions,  the  two  rotations  a and  8 are  related  to 
the  lateral  deflection  w by  the  derivatives 

flw 

rv  = — — 


0y 


R - 9W 
B ‘ ' 0X 


(14) 

(13) 


For  the  pure  bending  element,  the  matrices  u and  R become 

U . {»1o.161»jaj6Jwfeo|!6few1o,161}  (16) 

R - (17) 


9 


The  rectangular  bending  element  chosen  Is  described  In 


References  5 and  6 and  referred  to  in  the  latter  reference  as 

I 

the  Melosh  or  M element.  In  bending,  this  element  Is  assumed 
to  distort  along  the  edges  with  shapes  defined  by  beam  dis- 
placement functions.  These  displacements  are  assumed  to  de- 
crease linearly  toward  the  opposite  edge.  In  addition  it  is 
assumed  that  twisting  effects  are  resisted  by  a state  of 
uniform  twist  in  the  element.  This  state  distorts  the  element 
to  a simple  hy pe rb o 1 i c- p a r ab o 1 o i d shape.  Plots  documenting 
the  excellent  convergence  characteristics  of  this  element  are 
presented  in  Reference  6 as  well  as  in  the  following  section 
of  this  report.  These  latter  plots  were  constructed  from  new 
data  obtained  with  the  use  of  the  FINEL  program. 

The  Appendix  shows  the  stiffness  matrix  k for  this  rect- 
angular bending  element  as  well  as  the  matrices  fe  and  k (in 
symbolic  form)  for  the  combined  membrane  and  bending  element 
(oriented  with  respect  to  the  element  coordinate  system  as 

shown  in  Figure  2).  I 

■ 

THE  TRIANGULAR  COMBINED  MEMBRANE  AND  BENDING  ELEMENT 

As  with  the  rectangular  combined  element,  the  stiffness 
matrix  for  the  triangular  combined  membrane  and  bending  ele- 
ment is  obtained  by  combining  the  stiffness  matrices  of  a 
triangular  membrane  element  and  a triangular  bending  element. 

The  membrane  element  employed  is  the  well  known  "constant 
stress"  triangle  included  in  the  original  FINEL  program  as 
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Element  1 and  described  in  Reference  2.  The  assumed  membrane 


stress  distribution  for  this  element  is  illustrated  in  Fig- 
ure 3.  The  triangular  bending  element  chosen  to  provide  the 
bending  stiffness  matrix  is  described  in  the  following  para- 
graphs . 

For  the  triangular  combined  element,  oriented  as  shown 
in  Figure  3,  the  stiffness  matrices  k and  R are  defined  in 
terms  of  the  node  displacement  and  node  force  matrices  repre- 
sented by 


u 

= {uiv1wia1BiujVjwjaj6jUfeVfewfeafeBfe} 

(18a) 

(1) 

u 

= {ujwjajBjukvfewkakV 

(18b) 

(2) 

u 

= {u.v.w.a.g.v.} 

i i l i i j 

(18c) 

and 

R 

= {XYZ.K.M.X.YZ.K.MX/YiZ/KlMi} 
iiiiijjjjjfefefefefe 

(19a) 

R(1) 

= {X.Z . K , M . X , Y.Z.KiM,} 

(19b) 

r(2) 

= {X  Y.Z.K.M  Y.  } 
i l i i i j 

(19c) 

For 

the 

pure  bending  element,  the  u and  R matrices 

become 

(20) 

! ‘ (ZiKlMiVjVfcKtMt1 

(21) 

The 

particular  triangular  bending  element  selected  is  a 

nonconf  orming 

element  which  is  described  in  detail 

in  Rrfer- 

ences  3 

and  7 . 

The  displacement  function  for  this 

e 1 emen  t 

uses 

the 

so-called  "area  coordinates"  to  ensure  symmetry  and 

11 
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simplicity  of  statement.  The  formulation  of  this  function  is 
based  on  the  hypothetically  necessary  condition  (for  conver- 
gence) that  the  displacement  function  is  capable  of  represent- 
ing constant  curvature  (strain)  states  throughout  the  finite 
element,  irrespective  of  element  size  or  shape. 

The  actual  derivation  of  the  stiffness  matrix  k for  this 
element  is  carried  out  by  first  determining  a stiffness 

* 

matrix,  labeled  k , for  a simply  supported  (at  the  nodes)  ele- 
ment, and  then  transforming  this  matrix  to  the  desired  stiff- 
ness matrix.  This  transformation  is  carried  out  by  means  of  a 

| 

transformation  matrix  T which  is  a function  of  the  element 
geometry  and  its  orientation  with  respect  to  the  local  element 

I 

coordinate  system.  The  expression  relating  k , T,  and  k is 

k = TT  • k*  • T (22) 

* 

Because  of  the  complexity  of  the  components  of  the  k matrix, 

~ 

the  matrix  k is  computed  numerically  in  the  FINEL  program  by 
performing  the  matrix  multiplications  indicated  by  Equa- 
tion (22). 

* 

The  Appendix  contains  the  k and  T matrices  for  the  tri- 
angular bending  element  oriented  as  shown  in  Figure  3 as  well 
as  the  matrices  k and  k (in  symbolic  form)  for  the  triangular 
combined  element. 
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BENDING  CONVERGENCE  CHARACTERISTICS 
OF  THE  COMBINED  ELEMENTS 


In  order  to  evaluate  the  bending  behavior  of  the  two  com- 
bined elements,  a series  of  analyses  was  carried  out  on  a uni- 
form square  plate  laterally  loaded  by  a concentrated  force  at 
its  center.  Two  edge  support  conditions  were  considered:  all 

edges  simply  supported,  and  all  edges  clamped.  This  problem 
was  chosen  because  of  the  availability  of  the  "exact"  solu- 
tions in  a number  of  texts. 

In  this  study,  the  deflection  of  the  center  of  the  plate 
w^  was  taken  as  the  measure  of  the  quality  of  the  approxima- 
tion. This  deflection  is  represented  by  a nondimensional 
"deflection  coefficient,"  designated  C^,  which  is  defined  as 


Cd  = 


w D 

c 

Qs2 


(23) 


where  D is  the  flexural  rigidity  of  the  plate, 
s is  the  edge  dimension  of  the  plate,  and 
Q is  the  magnitude  of  the  concentrated  force. 

Using  linear  plate  bending  theory,  the  "exact"  solutions  for 
the  deflection  coefficient  C ^ for  the  two  cases  of  edge  sup- 
port are  found  to  be 

C ^ = 0.01160  Simply  Supported 


and 

C,  = 0.00560  Clamped 
d 
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Because  of  its  double  symmetry,  only  one-quarter  of  the 


plate  was  considered  in  the  finite  element  analysis.  The 
convergence  characteristics  of  the  two  elements  in  bending 
were  studied  by  using  a series  of  different  element  mesh 
arrangements  and  sizes  in  the  analysis  of  each  case.  The  mesh 
size  is  indicated  by  a mesh  number  r which  refers  to  the  num- 
ber of  rectangles  along  the  side  of  the  c omp 1 e t e plate.  The 
three  basic  mesh  arrangements  employed  in  this  study  are  shown 
in  Figure  4 for  a value  of  the  mesh  number  r = 8.  Figure  5 
shows  another  mesh  arrangement,  in  this  case  for  r = 4,  for 
which  two  additional  data  points  were  plotted.  This  mesh  was 
investigated  primarily  to  allow  for  additional  comparisons  of 
data  produced  using  the  FINEL  program  with  that  published  in 
Reference  7.  As  is  noted  in  that  reference,  however,  this 
mesh  arrangement  does  not  converge  to  the  correct  solution 
when  used  with  nonconforming  elements. 

The  results  of  this  series  of  analyses  are  presented  in 

Figures  6 and  7.  These  figures  are  plots  of  the  nondimen- 

sional  deflection  coefficient  C,  versus  the  mesh  size  r. 

d 

Both  Figures  6 and  7 illustrate  the  excellent  bending 
convergence  characteristics  of  the  rectangular  combined  ele- 
ment. Other  known  results  obtained  with  this  element  in  its 

pure  bending  form  also  demonstrate  excellent  convergence.^ 
Comparisons  with  these  older  data  suggest  that  the  drop  in  the 
curve  for  Mesh  A at  r = 16  on  Figure  6 is  the  result  of  pro- 

These  examples  indicate  that  the  use  of 


gram  round-off  error. 
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J 


this  rectangular  element  can  give  good  results  for  plate  bend- 
ing problems  even  with  relatively  coarse  meshes. 

It  is  also  apparent  from  the  figures  that  good  results 
for  plate  bending  probl  "ns  ar  possible  using  the  triangular 
combined  element,  but  the  .ale  of  convergence  is  very  much  a 
function  of  the  mesh  arrangement.  Discussions  of  convergence 
problems  with  nonconforming  triangular  elements  indicate  that 
certain  mesh  arrangements  do  not  converge  to  the  correct 
3 , 7 

solution.  ’ For  the  arrangements  considered  here,  these 
sources  indicate  that  Meshes  B and  C will  converge  to  the  cor- 
rect solution  whereas  Mesh  D will  not.  Readers  are  urged  to 
examine  the  references  cited  for  a more  detailed  discussion 
of  this  subject.  This  behavior  suggests  that  for  best  results 
with  plate  bending  problems,  rectangular  elements  should  be 
used  whenever  feasible.  Should  the  geometry  of  the  problem 
be  such  that  the  use  of  triangular  elements  is  either  desired 
or  unavoidable,  it  is  essential  that  the  user  be  aw a r e of  the 
convergence  problems  associated  with  their  use. 

SUMMARY 

Three  new  elements  have  been  added  to  the  library  of  the 
FINEL  finite  element  program:  a rectangular  membrane  element, 

a rectangular  combined  membrane  and  bending  element,  and  a 
triangular  combined  membrane  and  bending  element.  The  rect- 
angular membrane  element  is  a special  case  of  the  general 
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quadrilateral  membrane  element  included  in  the  original  FINEL 
program. 

The  rectangular  combined  element  exhibits  excellent  con- 
vergence characteristics  in  bending.  Good  results  for  plate 
bending  problems  can  be  achieved  with  this  element,  even  for 
relatively  coarse  meshes.  This  element  is  recommended  for 
such  problems  whenever  the  problem  geometry  permits. 

Good  results  for  plate  bending  problems  are  also  possible 
with  the  triangular  combined  element,  but  the  rate  of  conver- 
gence is  very  much  a function  of  mesh  arrangement.  Certain 
mesh  patterns  do  not  converge  to  the  correct  solution.  An 
understanding  of  the  convergence  problems  with  nonconforming 
bending  triangles  is  necessary  for  users  of  this  element. 


LOCAL  ELEMENT 
COORDINATE  SYSTEM 


z 


DISPIACEMENT  COMPONENTS 

Translational  - u,  v,  w 
Rotational  - OC , (3  , T 

LOAD  COMPONENTS 

Force  - X,  Y,  Z 
Moment  - K,  M,  N 

Figure  1 - Element  Coordinate  System 

J 
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Figure  2 - Rectangular  Element  Geometry  and  Assumed  Membrane 

Stress  Distribution 


Figure  3 - Triangular  Element  Geometry  and  Assumed  Membrane 

Stress  Distribution 
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MESH  B 

Triangular  Elements 


MESH  C 

Triangular  Elements 


Figure  5 - Special  Mesh  D (r  - 4) 
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DEFLECTION  COEFFICIENT 


Figure  6 - Deflection  Coefficient  versus  Mesh 
Size  for  Simply  Supported  Plate 
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DEFLECTION  COEFFICIENT 


APPENDIX 


ELEMENT  STIFFNESS  MATRICES 


RECTANGULAR  MEMBRANE  ELEMENT 


The  stiffness  matrix  fe  of  the  rectangular  membrane  ele- 
ment is  computed  from  a matrix  product  given  by 

T 


(B  l) 


B 


-1 


(ID 


-1 


For  this  rectangular  element,  Matrices  B,  B , and  J have  the 
following  forms: 


B = 


1 

E 


-vyfe 


~2XjXfe 


7V 


■vyfe  -1vyk 


-vx 


j 


yfe 


xj  2Xjyfe  ~Vxj 


yfe 


i 

Tyfe 


2 2 


2 ( 1 +v ) y 


-j(vx  . + y f,  ) 2 ( 1 + v ) y^ 

x.y,.  0 


If  the  matrix  multiplications  indicated  by  Equation  (11) 


are  performed,  the  stiffness  matrix  fe  is  the  result: 


fe 


( C1 i + C3  2 > 

(”C2l"C32) 

(C11+C32) 

SYMMETRIC 

|(l+v) 

“(1+V) 

(C12+C31> 

(C2l"C32) 

<-Cll+C32) 

d-3v) 

(C11+C32) 

— j( l-3v  ) 

|(l-3v) 

(C22_C31) 

f(l+V) 

where 


'11 


'12 


'21 


'22 


'31 


(yfe/Xj) [i+j(i-v2) ] 

(x  /y  ) [i+I(i-v2>] 

3 k 3 

(y,  /x.) [ l-4(l-v2 ) ] 

k 3 J 

(*j/yfe)[ l-j(l-v2) ] 

^yfe  ^xj  * ^ 1_V^ 
kx./yWl-v) 


Matrix  P is  required  in  order  to  obtain  the  stiffness  matrix 
k from  k as  indicated  by  Equation  (8).  This  matrix  is  given 


r-H 
1 ' 

0 

1 

0 

0 

0 

0 

0 

-1 

■Vxj 

0 

Vxj 

1 

0 

0 

0 

0 

-1 

0 

0 

0 

1 

0 

0 

-1 

-yfe/xj 

0 

Vxj 

0 

0 

1 

0 

O 

1 

0 

0 

-1 

0 

0 

0 

1_ 

Performing  the  operations  indicated  by  Equation  (8)  will  pro- 
duce the  matrix  k-  This  matrix  may  be  written  in  symbolic 
form  as 


kll  k 1 2 

k 13 

k 14 

k 15 

k 1 6 

k 1 7 

k 1 8 

k 2 2 

k 2 3 

k 24 

k 2 5 

k 2 6 

k 2 7 

k 28 

fell 

k 34 

fe  12 

fe13 

14 

fe15 

k 4 4 

k 4 5 

k 46 

k 4 7 

k 4 8 

fe22 

fe23 

k24 

k25 

SYMMETRIC 

k33 

j 4 

fe35 

b44 

k45 

k55 

in  which  it  may  be  noted  the  entire  k matrix  appears.  This 
demonstrates  how  the  k stiffness  matrix  can  be  determined 
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from  the  k stiffness  matrix  simply  by  inspection  when  u 
and  u(2)  are  known. 


RECTANGULAR  BENDING  ELEMENT 

The  stiffness  matrix  k for  the  rectangular  bending  ele- 
ment is  defined  in  terms  of  the  u and  R matrices  described 
earlier  as  Equations  (16)  and  (17).  This  matrix  can  be  more 
conveniently  written  as  the  sum  of  four  considerably  simpler 
matrices 


— — {k  +k,+k  + k , } 

ab  ~ a ~b  ~c  ~ d 


where  a = x . 

J 


yk 


D = 


Et 

12 ( 1-v  2 ) 
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I 7? 


6 

3b  2 
0 0 0 

0 0 0 6 


SYMMETRIC 


0 0 

0 

3b  2b2 

9 

0 0 

0 

0 0 

0 

L. 

a 

'-b 

b2 

-6  -3b 

0 

0 0 

0 

6 

3b  b2 

0 

0 0 

0 - 

3b  2b2 

0 0 

0 

0 0 

0 

0 0 0 

0 0 

0 

-6  -3b 

0 

0 0 0 

6 

0 0 

0 

3b  b2 

0 

0 0 0 - 

3b 

2b2 

_0  0 

0 

0 0 

0 

0 0 0 

0 

0 

" 72 

30b 

0 

-30a 

-25 ab  0 

SYMMETRIC 

-72 

-30b 

6a 

72 

-30b 

0 

5 ab 

30b 

0 

-6a 

- 5 ab 

0 

30a 

2 5 ab 

0 

-72 

-6b 

30a 

72 

6b 

6a 

72 

6b 

0 

- 5 ab  - 

6b 

0 

ab 

-30b  0 

30a 

5ab 

0 

6 a 

-ab 

0 

-30a  25ab 

0 

72 

6b 

-6a 

72 

-6b  - 

30a 

-72  30b 

6a 

72 

-6b 

0 

ab 

6b 

0 5ab 

30b  0 

5 ab 

-30b 

6a 

ab 

0 

30a 

-Sab 

0 

-6a  5ab 

0 

30a 

0 

-2  5 ab 
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kd  = 2 ( 1-v ) 


TRIANGULAR  BENDING  ELEMENT 

It  was  stated  earlier  that  the  stiffness  matrix  k of  the 
triangular  bending  element,  defined  in  terms  of  the  u and  R 
matrices  presented  as  Equations  (20)  and  (21),  is  computed 
from  the  matrix  product  given  by 

k = TT  • k*  • T (22) 
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P 


For  the  orientation  of  the  triangle  as  defined  in  Figure  3, 
the  transformation  matrix  T is  the  following 


T 


1 

2 A 


(d-e ) 
-f 

(d-e) 

-f 

(d-e) 

-f 


2 A 0 e 0 0 -d  0 0 
0 2A  f 0 0 00  0 

0 0 e 2A  0 -d  0 0 

0 0 f 0 2A  0 0 0 

0 0 e 0 0 -d  2A  0 

0 0 f 0 0 0 0 2A_ 


where  A = -iyx.y,  the  area  of  the  triangle, 
i 1 fe 

d = Xj 

e = x,  , and 
fe 

f = V 

* 

The  stiffness  matrix  k for  the  "pinned"  element  can  be  cal- 
culated by  following  the  procedures  described  in  References  3 
and  7.  The  general  term  of  this  matrix  can  be  most  conven- 
iently represented  as  a function  of  three  arbitrarily  named 
matrices:  F,  L,  and  S. 
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DI<Fll<P1J+vF2J)+F21(F2j+VFlJ) 


1] 


+7<1-V>F31F33  )A+{Fli  aij+VI'2j  ) + Lll  <FlJ+UF2J  > 


+F2i(L2j+VLlj)+L2i(F2j+VFlj)+2(1"V)(F3iL3j 

tL3iF3J),r;<d+'s)  + fFll<SlJ+VS2i)  + Sll(Fl) 
+uF2J>+F21<S2J+VSl1>+S21<F2j+VFl1> 


+l(1-v,<F3iS3j+S31F3j),1f  +(Lll<Llj+VL2J> 

+ L2i(L2j+VLli  )‘lT(1_V)L3iL3j  >ff(d2+de+e2) 
+{Lli(Slj+VS2j)+Sli(Llj+VL2j)+L2i(S2j+VSlj) 


+S21<L23+VLlJ>+7<1-v)<L31S3J+S31L3j)fl  <d+2e) 

+(Sli(Slj+VS2j)+S2i(S2J+VSlJ) 


+2 ( 1-v) S31S3j >12  ] 

The  matrices  F,  L,  and  S are  defined  as  follows 


F = 


8 A' 


2 3 
-4d  f 


3 2 
-2f  d 


,2,2 


,2  ,2 


2,2 


d f ( 4d  -d  ef(e  d ef  d ef(2d  2d  f (d  d 


■ 3 e ) 


+ d) 


-3e) 


,2,3  ,2,2 


2,3  ,2,2 


-3d  fJ  d f (3d  d f d f (2d  - 2 d “ f 


-5e) 


- 5e  ) 


-e) 


,2,3 
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0 

2ef ( 2e 
-3d) 

2 2 
f (2e 

-d) 


0 

6 d f 3 

0 

6df  3 

0 

0 

df 2 (2e 

df  (4e2 

d f 2 (2e  f d ( 4 e 2 

2 d f 2 ( 2 e 

df  (d2 

-3d) 

-3de+d2) 

+d)  - 

5 de  + 2 d 2 ) 

-d) 

+4de-4e2) 

2 d f 3 

2 d f 2 ( 5e 
-3d) 

2 f 3 d 

2 f 2 d ( 5 e 
-2d) 

4d  f 3 

2df 2 (d 
-2e)  _ 

-df3 

df 2 (3d 
- 5e  ) 

-f3d 

f 2d ( 2 d 
-5e) 

- 2 d f 3 

df 2 (2e 
-d) 

3d  f ( d 

3de (d2 

-3de  f (e 

3d  e ( e 

6df (de 

3de ( e- 

-e)(e-2d)  -e2) 

+d) 

-2d) (d- 

, A2  2 

e ) -d  -e 

) d) ( 2e-d) 

2df 2 (3d 

2df ( 3de 

-2df  2 (2 

e 2fd(5de 

4df 2 (d 

2 d f ( 4 e 

-2e) 

2 2 
-4e  -d 

) +d ) - 4 e 2 - 2 d 

2)  - 2e  ) 

- 4 d e - d ) 

Because  of  the  complicated  algebraic  form  of  the  com- 

"k 

ponents  of  k , it  is  necessary  in  the  FINEL  program  to  compute 
the  stiffness  matrix  k for  this  element  numerically  within  the 
program  according  to  the  matrix  operations  dictated  by  Equa- 
tion (22). 


COMBINED  MEMBRANE  AND  BENDING  ELEMENTS 

As  discussed  earlier,  the  stiffness  matrices  for  the  com- 
bined membrane  and  bending  elements  are  formed  by  combining 
the  stiffness  matrices  of  the  appropriate  membrane  and  bending 
elements.  The  following  notation  is  employed  to  distinguish 
in  the  combined  stiffness  matrices  between  stiffness  terms 
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associated  with  the  membrane  stiffness  matrix  and  terms  asso- 


ciated with  the  bending  stiffness  matrix: 

mlj  ^ij  ^membrane 

b = (k  ) 
ij  VKij  'bending 

Rectangular  Combined  Membrane  and  Bending  Element 

Using  the  notation  just  defined,  the  stiffness  matrix 
for  the  rectangular  combined  membrane  and  bending  element, 
defined  in  terms  of  the  u and  R matrices  presented  as  Equa- 
tions (12a)  and  (13a),  becomes 
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I 7? 
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10,10 


I FT 


By  inspection,  the  k stiffness  matrix  for  this  element  is 


'54  “55 
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0 
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Triangular  Combined  Membrane  and  Bending  Element 

Here  the  u and  R matrices  are  defined  by  Equations  (18a) 
and  (19a),  respectively. 
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By  inspection,  the  corresponding  k stiffness  matrix  for 
the  triangular  combined  membrane  and  bending  element  is  given 
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